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O ■ ABSTRACT 
<N 

Context. We present a detailed spectroscopic analysis of the stars with helium-dominated spectra in the ESO Supernova la Progenitor 
Survey (SPY). 

Aims. Atmospheric parameters, masses, and abundances of trace hydrogen are determined and discussed in the context of spectral 
evolution of white dwarfs. 

^ , Methods. The spectra are compared with theoretical model atmospheres using a fitting technique, leading to determinations of 

5^ . effective temperature, surface gravity, and hydrogen abundance. 

' Results. Our final sample contains 71 objects, of which 6 are new detections and 14 are reclassified from DB to DBA because of the 

' presence of H lines. One is a cool DO with weak He ii lines, 2 are composite DB+dM. 55% of the DB sample show hydrogen and are 

thus DBA, a significantly higher fraction than found before. 

Conclusions. The large incidence of DBA, and the derived total hydrogen masses are compatible with the scenario that DBs "reap- 
pear" around 30 000 K from the DB gap by mixing and diluting a thin hydrogen layer of the order of 10 '"^ Mq. This hydrogen mass 
^ ' is then during the evolution continuously increased by interstellar accretion. There are indications that the accretion rate increases 

, smoothly with age or decreasing temperature, a trend which continuous even below the current low temperature limit (Dufour 2006). 

CO . A remaining mystery is the low accretion rate of H compared to that of Ca observed in the DBZA, but a stellar wind extending down 

0^ ' to the lowest temperatures with decreasing strength might be part of the solution. 

O 

. 

1 . Introduction While the properties of the hydrogen-atmosphere DA white 

• • dwarfs as a group have been closely investigated in numerous 

_ > About 20% of all white dwarfs (WD) have atmospheres in which studies, the DBs received less attention. They are less numerous 

■k> helium is the main constituent. Such hydrogen-deficient stars are than the DAs and thus it is more difficult to assemble data sets 

rN believed to form when a cooling stellar remnant experiences that allow meaningful statistical analyses. The most comprehen- 

Cj a very late helium-shell flash and goes through a born-again give study of DB stars so far was that of a sample of 80 DBs, 

episode, i.e., a last thermal pulse, in which almost all residual carried out by Beauchamp et al (1996 ) 

hydrogen is burnt (see, e.g., Althaus et al. gOOSl and references ^.^^^ Beauchamp's study, the number of known DB stars 
^' was not significantly increased until two surveys, the Sloan 
As these stars advance on the WD cooling sequence, they Digital Sky Survey (SDSS, York et al. 2000) an d the Supernova 
can initially be observed at temperatures down to approxi- ja Progenitor Survey (SPY, Napiwotzki et al. |2003l l collected 
mately 45 000 K as hot DO stars with spectra dominated by He ii new samples of DB spectra. The SPY sample contains high- 
fines, and subsequently as DB stars whose spectra show only resolution spectra of about 60 DB white dwarfs. These data al- 
Hei, at temperatures that are mostly below 30 000 K. Finally, Iqw to determine well-constrained parameters of a large num- 
when Helium becomes spectroscopicaUy invisible below effec- ber of DB stars for first time since Beauchamp's work. The 
tive temperatures around 10 000 K, the helium atmospheres dis- sDSS spectra, on the other hand, are of medium resolution, but 
play either a featureless DC-type spectrum, or a spectrum dom- ^uch more numerous. Thus, the SDSS detected a large num- 
inated by the features of trace constituents, i.e., dredged-up car- ber of new DBs. The latest version of the catalog of WDs in the 
bon leading to DQ spectra or accreted metals leading to a DZ- sDSS (Eisenstein et al. 2006a), based on the SDSS data release 
type spectrum. Besides these main types, DB spectra which 4 (dR4), contains 713 DB stars. About 560 of these stars are or- 
show traces of hydrogen or metals are known, classified as DBA dinary DBs which only exhibit He i in their spectra. The second 
or DBZ stars. largest subgroup is that of DBA stars with weak Balmer lines 
in the otherwise He i-dominated spectrum; about 90 such stars 
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* Based on data obtained at the Paranal Observatory of the European are found in the SDSS sample. Many of these newly detected 
Southern Observatory for programmes 165.H-0588 and 167.D-0407. DBs are faint, which makes detailed spectroscopic studies diffi- 
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cult. Nevertheless, the almost tenfold increase of the number of 
DB stars has allowed the first statistically meaningful analysis of 
some key questions connected to the evolution of DB stars: 

As a helium-atmosphere WD cools, it encounters two tem- 
perature regions of special interest, the "DB gap" between tem- 
peratures of 45 000 K and 30 000 K, and later the instabihty re- 
gion of the V777 Her variables near temperatures of 25 000 K. 
The DB gap region has been enigmatic for some time since no 
DB stars were found in this range of temperatures, raising ques- 
tions about the very nature of helium-atmosphere white dwarfs. 
However before the SDSS, only few cool DO and hot DB stars 
that constrained the gap were known and thus the reality of the 
gap was not well established. Eisenstein et al. (12006b') found 
several hot DB stars from the SDSS within the DB gap, some of 
which exhibit weak He n lines and are thus situated in the higher- 
temperature part of the DB gap. The gap has however not been 
closed by their discovery because they also find that the number 
density of DB stars in the gap is still lower by a factor of about 
2.5 than what would be expected from the density of cooler DB 
stars. 

As a possible explanation, it has been speculated that cool- 
ing DA stars might change into DBs, provided that their hy- 
drogen layer is very thin, on the order of Mu ~ 10^'^Mq. In 
this case, the superficial hydrogen would be convectively mixed 
into the more massive Helium layer when the atmosphere cools 
to temperatures lower than 30 000 K and the Helium layer be- 
comes convective (see Eisenstein et al. I2006bl and references 
therein). Some peculiar objects of mixed helium/hydrogen at- 
mosphere compositions, the DAB and some of the DBA stars, 
might be objects which are currently undergoing such a transi- 
tion and could therefore provide insight into the viability of this 
scenario. 

The majority of cool DBA stars however can more easily be 
explained as DB stars that are accreting hydrogen from the ISM. 
They were long found to comprise about 20% of the DB stars, a 
fraction that was first estimated by Shipman, Liebert & Green 
( 1 1987 J and later confirmed, e.g, by Beauchamp et al. (19961). 
Similarly, 22% of all DB stars in the catalogue of McCook & 
Sion ( fT999b are classified as DBA, and in the recent SDSS DR4 
sample, the fraction is about 13%. It is clear that the fraction 
of DBA among the DB stars is a question of the observational 
hydrogen detection limits. 

It has long been expected by some authors that in fact all 
DBs accrete small traces of hydrogen and would thus appear as 
DBA if only small enough traces of hydrogen could be detected. 
Most studies so far, including the SDSS, could only detect hy- 
drogen abundances, relative to helium, larger than a few times 
10"^. These abundance measurements were mostly based on the 
appearance of Hy and Hfi in the spectra. Lower abundances can 
be identified from the Lya line in UV data (e.g. Provencal et al. 
I2OOOI 1. but this method is more involved and not easily applied to 
large samples. Hunter et al. ( 120011 ) demonstrated that also opti- 
cal spectra allow a measurement of lower hydrogen abundances 
if they cover the Ha line. In their sample of 24 DB objects, the 
DBA fraction is raised to about 30%. With observations of a bet- 
ter S IN and higher spectral resolution, this fraction was expected 
to be even larger, with the expectation that probably hydrogen 
might be identified in a majority of the DB stars. 

Fontaine et al. (120051) proposed that stellar winds of WD stars 
might prevent hydrogen accretion at temperatures in excess of 
20 000 K, thus enabling the conversion from DB to DBA stars 
only at lower temperatures. Observations with increased sensi- 
tivity that allow to detect increasingly lower H masses in DB 



stars should allow to test this effect and to provide new insight 
into the evolution of DBA stars. 

The high resolution WD spectra that were obtained by SPY 
contain the largest sample to date of high-quality DB spectra and 
it should thus be possible to determine lower hydrogen abun- 
dances than was possible before from optical data. 

The SPY data will also be helpful to investigate the ques- 
tion if invisibly small amounts of trace hydrogen are present in 
hot DB stars with temperatures above 20 000 K. This question is 
of importance because even such small traces of hydrogen sig- 
nificantly affect the shape of the He i lines and thus complicate 
the determination of temperatures from the spectra. Beauchamp 
et al. (1999,) compared temperature estimates based on pure He 
atmospheres to estimates including effects of invisibly small hy- 
drogen contaminations, and found a shift of up to 3000 K in tem- 
perature. The temperatures of hot DBs are therefore not well 
constrained so far. An improvement of these temperature esti- 
mations would greatly benefit the analysis of the V777 Her in- 
stability strip. 

2. The SPY data 

SPY is a radial velocity survey that was conducted to test the 
double degenerate channel of the formation of supernovae la. 
About 800 white dwarfs were observed in the course of the sur- 
vey, assembling a large collection of high-quality white dwarf 
spectra. Most of the DB stars among the SPY objects were se- 
lected as targets based on the WD catalog (McCook & Sion 
'19991) and from a study of DBs in the Hamburg ESO Survey 
by Christlieb et al. ( I1997I ). A more detailed discussion of the tar- 
get selection for the SPY observations is given in Napiwotzki et 
al. (2001,2007). 

The SPY spectra were employed for a number of studies be- 
yond the original scope of SPY, of which the work of Koester et 
al. ( 120011) , hereafter "paper I", is most relevant to our analysis. 
They have derived temperatures and gravities from a preliminary 
sample of about 200 SPY stars, including 13 DB stars. 

Since they have already given a detailed description of the 
data properties, we only briefly summarize the observations and 
data reduction: 

The spectra were obtained with UVES, a high resolution 
echelle spectrograph at the ESO VLT telescope. UVES was used 
in a dichroic mode, resulting in small gaps, k %QA wide, at 
4580 A and 5640 A in the final merged spectrum. The spectral 
resolution at Ha is R = 18 500 or better, and the S /N per binned 
pixel (0.05 A) is S/N^ 15 or higher. 

The spectra were reduced with the ESO pipeline for UVES, 
including the merging of the echelle orders and the wavelength 
calibration. Paper I found that the quality of these automatically 
extracted spectra is very good, except for a quasi-periodic wave- 
like pattern that occurs in some of the spectra. This has since 
been largely removed by additional processing by collaborators 
of the SPY project at the University of Erlangen-Niirnberg. A de- 
tailed description of the reduction procedure applied to the spec- 
tra analyzed here will be published in Napiwotzki et al. (20071 ). 
Some artifacts remain in the data, but they do not significantly 
affect the spectral analysis. 

For our analysis, the data have been rebinned to a step size 
of 0. 1 A. This is in contrast to the analysis of paper I which used 
a binning of 1 A. Therefore, we include the paper I subsample in 
our analysis, now with the smaller step size. Further details about 
the observational setup and the data reduction can be found in 
Napiwotzki et al. (2003 ) and in paper I. 
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3. Analysis method 

3.1. Model atmosphere fits 

We fit the data with model spectra from six different grids of 
helium atmosphere models with different hydrogen abundances, 
from H/He = 10"^ down to H/He = 10"^, and with a grid of 
pure He atmosphere models. The grids cover a range of temper- 
atures from lOOOOK to 50 000 K and logarithms of the surface 
gravities from 7.0 to 9.0. The input physics of the models are de- 
scribed in Finley et al. ( 1997), Homeier et al. ( 1998), Koester & 
Wolff (2000 ), and Koester et al. ( 2005a ). We use ax^ minimiza- 
tion fitting routine that is based on the Levenberg-Marquardt al- 
gorithm (Press et al.'1992^ to derive the best fitting eff'ective tem- 
perature and surface gravity for each spectrum. Some more de- 
tails on the fitting process can be found in Homeier et al. (|1998l l. 

The SPY spectra are not reliably flux calibrated, and there- 
fore we do not use their continuum slopes for the fitting but only 
the profiles of the helium lines. We determine a reference con- 
tinuum level at each wavelength to which the models are scaled; 
this reference level is derived by interpolating the flux levels 
measured next to the He i lines. 

The Helium lines reach their maximum strength around 
Teff = 20 000 K, and therefore two solutions for the fit of a spec- 
trum are normally possible, one on either side of the line strength 
maximum. We employed two T^ff starting values for the fitting, 

15 000 K and 25 000 K, to cover all possible solutions. 

For most objects it is not difficult to pick the correct solu- 
tion; if the temperature difference of both solutions is large, the 
line profiles are different enough to allow an easy decision for 
the best fit. Furthermore, for many objects there is only one so- 
lution, either because both fits converge to a single solution, or 
because the fit from the hot starting value converges to exceed- 
ingly high values of the temperature. For the remaining few ob- 
jects, a careful study of the fit quality was necessary to determine 
which solution is to be preferred; for some objects the hydro- 
gen abundance or the presence of pulsations provided additional 
constraints. 

There are two stars in our sample that show composite spec- 
tra. The red part of the spectra was not used in the fits of these 
DB-t-dM binaries; only the Hei lines at short wavelengths of 
A < 4500 A were fitted, where the flux of the late-type com- 
panion is negligible. 

Below effective temperatures of approximately 15 500 to 

16 500 K - depending on surface gravity - neutral He van der 
Waals broadening becomes more and more important for the 
line shapes. Since our original model grids included only Stark 
broadening, we have calculated a second grid up to 18 000K 
which includes Stark broadening in a simplified form and van 
der Waals broadening combined. This greatly improved the 
agreement with the observed line profiles and strengths. Because 
the fit was not perfect, and also because in general the depen- 
dence of line profiles on surface gravity is smaller for the cooler 
DBs, we have decided to keep \ogg fixed at 8.0, in order to get a 
more robust estimate of the effective temperature. A fit in which 
log g is allowed to vary would converge to very high gravity val- 
ues. 

There are six objects in our sample whose spectra only show 
a strongly broadened 5876 A line, or only this line together 
with about equally strong Ha. These are very cool objects with 
helium-rich atmospheres; they are related to peculiar stars like, 
e.g., HS 0146+1847 (Koester et al. 2005a), which are dominated 
by helium but appear hydrogen-rich because of their very low 
temperature. Since only one or two lines are available in these 



spectra, no meaningful fit would be possible and thus we es- 
timate the temperature and the hydrogen abundance from the 
equivalent widths of He i 5876 A, and, if present, Ha, assuming 
a fixed \ogg - 8.0. 

3.2. Determination of H/He 

We obtain seven different fit results for each observed spectrum, 
one fit result per assumed hydrogen abundance. If a visual in- 
spection of the observed spectrum does not reveal the presence 
of Ha, we adopt the pure helium atmosphere solution. 

Otherwise, if Ho' is unambiguously present, as shown in 
Figs. lA.ll through IA.3I we measure the equivalent width of 
the line in the observed spectrum as well as in each of the 
six model spectra. Then, we interpolate the model equiva- 
lent widths EW, model temperatures T^ff, and model gravi- 
ties \ogg as functions of the assumed hydrogen abundance 
H/He: £W(H/He), T^giH/He), and log§(H/He). With the ob- 
served equivalent width EWobs, we can determine H/Hcobs from 
the inversion of EW(HfHe). We then adopt reff(H/Heob.s) and 
log ^(H/Hcobs) as the temperature and gravity of the object. 

For most stars, the Ha equivalent width does not rise 
monotonously over the whole range of modeled abundances. 
Instead, it mostly rises monotonously up to H/He= 10"^ and 
decreases again for H/He= 10"^. The reason for this effect is 
that a fit of the H/He = 10"^ models to the data results in a high 
temperature, for which the equivalent width is small in spite of 
the high abundance. A given equivalent width thus often allows 
two solutions for the H abundance. For all but three of our ob- 
jects the low-abundance solution proves to be the better fit of the 
data. For the objects for which instead the high abundance so- 
lution is preferred, the fit quality difference of both solutions is 
small. Therefore it has to be noted that for these objects a second 
solution of almost identical quality exists. 

3.3. Stellar masses and hydrogen mass fractions 

The masses were determined from T^ff and logg, and the evolu- 
tionary mass-radius relations of Wood (1 19951 1 for "thin" layers, 
appropriate for DB white dwarfs. 

The depth of the convection zones in DB white dwarfs 
for given mass and effective temperature were calculated with 
our stellar structure code. This code was originally developed 
for the calculation of finite-temperature mass-radius relations 
and equilibrium models for pulsation studies (Koester 119781 
Dziembowski & Koester [l981| l. This code is not fully evolution- 
ary, but solves the static stellar structure equations; The energy 
equation is replaced by setting l(r) oc m(r). This is a sufficiently 
accurate approximation, since in the outer layers it results in 
I K L- const and in the interior the temperature becomes nearly 
constant because of the high heat conduction of the degenerate 
electrons. The input physics - equation of state and opacities - 
have been updated and the outer boundary values are determined 
from accurate DB atmospheric models. The parameterization of 
the mixing length version is ML2, a - 0.6 (see e.g. Fontaine et 
al. [T98T] Tassoul et al. 1990 or Koester et al. [T994l for an expla- 
nation of nomenclature). The mass-radius relation predicted by 
these models is in reasonable agreement with more sophisticated 
calculations (e.g. Wood [T995] l. 

This code predicts the total mass in the convection zone as 
a function of the mass and effective temperature of the star. By 
multiplying this with the hydrogen abundance we obtain the total 
hydrogen mass within the convection zone in solar masses. 
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Fig. 1. The relative differences of the temperatures that we derive 
for two spectra of the same object. 
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Fig. 2. The same as Fig. [T] for log g. 



4. Results 



The properties and fit results of the DB stars without visible hy- 
drogen in their atmosphere are listed in Table[T] The data for the 
DBA stars are shown in Table |2] 



4.1. Uncertainties 

The formal uncertainties from the fitting routine are mostly 
low, and normally amount to only a few times lOK. However, 
they do not include systematic errors, which might be introduced 
by the data reduction, or by the lack of a flux calibration of the 
data. Thus the true uncertainties of the fit results are generally 
underestimated by the formal uncertainties. A better estimate 
can be obtained by comparing the fit results of the two indepen- 
dent SPY spectra of each object, an approach that was earlier 
used by Liebert et al. 2005 in their study of DA white dwarfs. 
We have done this for 24 of our objects for which both spec- 
tra are not too different in quality and aiTive at mean differences 
of < AFeff / Teff >= 3.22 % and < A log § >= 0.066 . The large 
<Areff > is dominated by two outliers at high temperatures. If 
these are omitted, the results are < AT^ff / Teff >= 2.03 % and 
<Alogg>= 0.058 . These results are plotted in Figs. [T] and |2] 



Fig. 3. The differences of the temperatures that we and other 
authors derive for common objects, plotted against our temper- 
ature. The -symbols show the comparison to the results of 
Beauchamp et al. (1999 ), the 'x' -symbols that to Castanheira et 
al. (2006 ) and the asterisks plot the comparison to the values of 
Friedrich et al. (i2000i) . 
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Fig. 4. The same as fig. |3] for log g. 



For six of the DBA objects we have conducted the full analy- 
sis for both spectra and can thus also directly estimate the uncer- 
tainties of H/He and Mh; five of the six objects show very low 
differences of abundance and hydrogen mass, <Alog(H/He)>= 
0.02 and <Alog(MH)>= 0.025, i.e., the total H masses differ 
by, on average, 6%. These estimates are of course based on a 
small number of objects and thus only allow to judge the order 
of magnitude of the uncertainties. 

A comparison to the parameters that other authors derive for 
common objects looks less satisfying; the differences in temper- 
ature between our results and those of Beauchamp et al. dl 9991 1, 
Friedrich et al. (2000), and Castanheira et al. (|2006l l, shown 
in Fig. |3] are scattered by more than +10% and the results of 
Beauchamp seem systematically hotter by a few %, or about 
1500 K. The surface gravity differences, shown in Fig.H) scatter 
by over +0.15. It has to be noted that Castanheira et al.'s results 
are derived from UV data, and that they use only pure helium 
models. Therefore we only compare our results for DB with- 
out hydrogen to their data. Beauchamp et al. assumed hydrogen 
abundances, whereas we can use observed H lines. They assume 
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H/He= 10"^ for their DBA stars, and certain lower hydrogen 
abundances corresponding to their detection Hmits for stars that 
they do not identify as DBA. The systematic differences of our 
results to theirs may in part arise because for all of our DBA 
stars that are in common with Beauchamp's sample, we derive 
lower abundances than he assumed. 

4.2. Hot DB stars 

Low amounts of hydrogen are spectroscopically invisible in hot 
atmospheres but can influence the shape of the helium lines 
through a change of the stratification and line broadening and 
thus affect the result of a model atmosphere fit. To test this 
eff'ect we have fitted the DB stars wifli T^ff > 20 000 K also 
with models of atmospheres with an admixture of an invisi- 
bly small amount of hydrogen, H/He = 10"^. These fit results 
are also listed in Table [T] The surface gravities do not differ 
significantly to those derived from hydrogen-free atmosphere 
models; the temperatures are on average larger if hydrogen is 
mixed into the atmosphere. The average temperature difference 
is Areff = 23 1 K, and thus on the order of the uncertainties of 
the measurements. The unknown hydrogen abundance therefore 
adds some additional uncertainty to the atmosphere parameters 
of hot DB stars, but this effect is much less pronounced in our 
results than in those of Beauchamp et al. ( 1999 ) and Castanheira 
et al. ( 120061 ). This is not surprising because the higher resolu- 
tion of our spectra allows us to use an upper limit of the abun- 
dance of spectroscopically invisible that is much lower than 
that of the earlier investigations, who used upper limits between 
H/He = 10"^-^ (Castanheira et al.) and H/He = 10 (Beauchamp 
et al.). 

Our sample includes only two V777 Her variables, and thus 
our results do not allow any meaningful statement on the struc- 
ture of the V777 Her instabihty strip. 

4.3. Peculiar objects 

4.3.1. The Cool DO Star WD 21 42-1 69 

The spectrum of WD 2142-169 is dominated by He i lines, but 
also shows a weak He ii 4686 A line, as well as Ha, li/3, and Ca ii 
lines. In a classification determined by the relative strength of the 
spectral features this star would be a DBO; however in accor- 
dance with Eisenstein et al. (2006b ) we classify it as a DO rather 
than a DB or DBO. The presence of He ii indicates that the tem- 
perature can not be much smaller than 40 000 K, but the weak- 
ness of the line suggests that the temperature of WD 2142-169 
is lower than that of most other DO stars, almost all of which 
haveTeff >45 000 K. 

WD 2142 -169 is listed as a DB star of Jeff = 15 900 K in 
Koester et al. ('2005W, and in fact such a fit can be obtained at 
a moderate fit quality, and with a low hydrogen abundance of 
H/He = 10"^, but of course such a fit does not reproduce the 
Hen line. The best fit is achieved with a much higher hydro- 
gen abundance, H/He = 10"^, and results in a temperature of 
7"etf = 39 100 K and a gravity of \ogg - 7.74. This fit is shown 
in Figs. |5] and |6] Note that NLTE eff'ects, which are important 
in hotter DOs, disappear below 40 000 K (Dreizler & Werner 
11996) . The He ii line is slightly stronger in the fit than in the spec- 
trum, indicating that the atmosphere is perhaps slightly cooler 
than the fit result. The B aimer lines are not well reproduced by 
this fit; they are stronger in the observed spectrum than in the 
fitted model, and thus the H abundance is probably higher than 
fitted, H/He > 10"^. A fit with a model of a low hydrogen abun- 
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Fig. 5. The fitted part of the spectrum of WD 2142- 169, shown 
as a black line, and the fit in red, shifted upwards by a constant 
amount for clarity. 
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Fig. 6. The part of the spectrum of WD2142-169 that contains 
the He ii 4686 A and H/3 lines. The spectrum is shown as a black 
line, and the fit in red, shifted upward. The feature at 4690 A is 
caused by a defect of the CCD that appears in many SPY spectra. 

dance usually overestimates the temperature. Thus, both the He ii 
line and the hydrogen features point to a true temperature that is 
even slightly lower than 39 000 K. 

WD 2142-169 is therefore one of the coolest DO stars 
that are known and one of very few helium-rich objects that 
populate the DB gap. Only one similar object is known so 
far, SDSSJ074538. 1+3 12205 (Eisenstein et al. 2006b). Like 
WD 2142-169, it shows weak Hen and Balmer lines in a spec- 
trum that is dominated by Hei. Eisenstein et al. derive a tem- 
perature of 39 800+2000 K from the SDSS photometry, however 
they do not obtain a reliable fit of the objects' spectrum. 

4.3.2. WD 0453-295 

This apparent DAB white dwarf was studied by Wesemael et al. 
(1994) and found to be a DA+DB binary, since no single star 
atmosphere model fit is able to reproduce the flux distribution of 
this object. 

The SPY spectra of WD 0453-295 (see Fig. Q show vari- 
able line core shapes, due to the orbital motion of the compo- 
nents. Napiwotzki et al. ( 12005) found that the core of Hff is split 



Table 1. Properties and parameters of the DB stars without visible hydrogen. Two lines of data are given for the hot objects (T^ff > 20 000 K); the second line contains the result 
of a fit with H/He=10"^ models. 



Object 


RA 


DE 


mag 




Alias 


T,B [K] 




M. [Mo] 


Comments 


WDOl 19-004 


01:21:48.27 


-00:10:54.4 


V = 


16.0 


G 27 1-047 A, LP 587-44, SDSS 


16542 


8.0 




4 


WD0127-3I1 


01:29:56.11 


-30:55:09.6 


V = 


14.51 


HE 01 27-3 no, GD 1363 


11002 


8.0 




only He I 5 876 A^ 


MCT0149-2518 


01:51:59.62 


-25:03:15.1 


B = 


15.91 


PHL1201 


16 834 


7.940 


0.557 




WD 0158-160 


02:00:56.69 


-15:46:09.7 


B = 


14.38 


G272-B2B, MCT 


25518 


7.875 


0.542 
















25 259 


7.891 




1 


WD 0203-181 


02:05:23.98 


-17:53:25.8 


V - 


16.0 


HE 0203-1807, G 272-152 


10757 


8.0 




only He I 5876 A^ 


WD 0249-052 


02:52:15.55 


-05:02:32.5 


B = 


16.1 


HE 0245-05 14, KUV 


17 908 


7.992 


0.589 




HE 0308-5635 


03:09:47.84 


-56:23:49.6 


B = 


14.02 


WD 0308-565, BPM 17088 


22 849 


8.060 


0.639 
















22 822 


8.081 




1 


WD 0349+015 


03:51:51.65 


+01:39:47.6 


B = 


16.3 


KUV 03493+0131 


18 741 


7.889 


0.535 




HE 0417-5357 


04:19:10.04 


-53:50:45.7 


B = 


15.12 


BPM 17731 


18733 


7.942 


0.563 




HE 0420-4748 


04:22:11.38 


-47:41:42.3 


B = 


14.7 




27 288 


7.808 


0.512 
















27 245 


7.809 




1 


HE 0423-1434 


04:25:51.72 


-14:27:52.0 


B = 


16.21 




16904 


7.794 


0.481 


new object 


HE 0429-1651 


04:32:13.81 


-16:45:08.5 


B = 


15.82 




15 659 


7.659 




DB+dM; new object 


WD 0615-591 


06:16:14.48 


-59:12:28.1 


V = 


14.09 


L 182-61, BPM 18164 


16714 


8.017 


0.600 




WD 0845-188 


08:47:29.52 


-18:59:50.7 


V = 


15.55 


LP 786-06, L 0748-70 


17 566 


7.969 


0.575 




WD 0900+142 


09:03:31.31 


+ 14:00:49.2 


B = 


16.48 


PG 0900+ 142 


15 678 


8.0 




4 


WD 1004-178 


10:07:07.91 


-18:05:25.8 


V = 


15.6 


EC 10047-1750 


16357 


7.700 




DB+dM 


WD 1046-017 


10:48:32.65 


-02:01:11.3 


V = 


15.81 


GD 124, GR387 


15138 


8.0 




4 


WD 1144-084 


11:46:54.07 


-08:45:48.5 


B = 


15.95 


PG 1 144-085 


16468 


8.0 




4 


WD 1252-289 


12:54:55.30 


-29:11:54.6 


V = 


15.85 


EC 12522-2855 


21029 


7.97 


0.595 




WD 1326-037 


13:29:16.37 


-03:58:51.8 


V = 


15.6 


PG 1326-037 


22645 


7.999 


0.603 


2 














22192 


7.995 




1 


WD 1336+123 


13:39:13.63 


+ 12:08:29.6 


B = 


13.9 


LP 498-26 


16779 


8.0 




4 


WD 1428-125 


14:31:39.64 


-12:48:55.5 


V = 


15.98 


EC 14289-1235 


21 586 


8.167 


0.701 
















21657 


8.179 




1 


WD 1445+152 


14:48:14.37 


+ 15:04:49.6 


B = 


15.55 


PG 1445+153 


23 234 


8.012 


0.611 
















24063 


7.968 




1 


WD 1542+182 


15:44:19.46 


+ 18:06:43.6 


V = 


14.72 


GDI 90, EG 193 


23186 


7.966 


0.585 
















22902 


7.965 




1 


WD 1612-111 


16:15:23.96 


-11:18:29.5 


V = 


15.53 


GD 198, EG 194 


22143 


8.073 


0.645 
















22 338 


8.076 




1 


WD 1654+160 


16:56:57.56 


+ 15:56:26.6 


B = 


16.15 


PG 1654+160, V 824 Her 


25 447 


7.846 


0.527 


V777 Her variable 














25 501 


7.839 




1 


WD 2129+000 


21:32:16.39 


+00:15:13.3 


V = 


14.73 


LP 638-4, PHL 028 


14414 


8.0 




4 


WD 2144-079 


21:47:37.34 


-07:44:13.1 


V = 


14.82 


G 026-03 1,LTT 8702 


16518 


7.895 


0.532 


DBZ 


WD 2234+064 


22:36:42.06 


+06:40:16.4 


.V = 


16.26 


PG 2234+064, GR907 


20 890 


8.035 


0.620 
















19 864 


8.028 




1 


WD 2354+159 


23:56:34.74 


+ 16:15:39.8 


B = 


15.78 


PG 2354+159 


24428 


8.180 


0.713 


DBZ 














22 822 


8.218 




1 


WD 2354-305 


23:56:37.07 


-30:16:25.4 


B = 


16.26 




18 006 


7.883 


0.530 





' Fit with H/He = 10 ^ models to check the influence of possible admixtures of invisibly small amounts of hydrogen. 

^ Beauchamp et al. ( 119991 1 list WD 1326-037 as a DBA; however no H is visible in the SPY spectra. According to Bergeron (2007, priv. comm.), the Beauchamp data show no hydrogen features 
and thus this object appears as a DBA in their paper due to a mis-classification. 

^ Temperatures were determined from the equivalent widths of He i 5876 A. The surface gravity was assumed to be log g = 8.0. 
Fitted with the alternate set of models with van der Waals line broadening. The surface gravity was kept fixed at logg = 8.0 for these fits. 



Table 2. Properties and parameters of the DB stars with hydrogen features in their spectra. For some objects, two Hnes of data are given, for the analysis results of each of the 
two independent SPY spectra of the star. 



Object 


RA 


DE 




mag 


Alia.s 


log (H/He) 


TciT [K] 


logg 


M. [Mo] 


Age [10" yi-j 


Mh [M 


d] 


Comments 




WD 0000- 1 70 


00:03:31.64 


-16:43:58.6 


V = 


14.67 


G 266-032, L0793-018 


^5.61 


13 354 


8.0 


0.583 


3.179 


4.410x10" 




6 




HE 0025-0317 


00:27:41.71 


-03:00:58.3 


B = 


15.69 




-3.13 


18 570 


8.293 


0.777 


2.017 


4.765x10" 


12 






HEOl 10-5630 


01:12:21.16 


-56:14:27.3 


B = 


15.86 




-4.18 


18 483 


8.120 


0.666 


1.350 


5.535x10- 


13 


DBAZ 
















-4.19 


18465 


8.132 


0.673 


1.393 


5.585x10- 


13 






WD0125-236 


01:27:44.57 


-23:24:49.3 


y - 


15.42 


G 274-039, MCT 


-5.01 


16757 


8.065 


0.629 


1.780 


1.287x10- 


12 


1 




HE 02 15-0225 


02:17:32.67 


-02:11:15.0 


B = 


16.13 


(PB6822) 


-5.67 


16776 


7.966 


0.571 


1.402 


3.383x10- 


13 


new object 




WD 0300-013 


03:02:53.20 


-01:08:34.6 


V = 


15.56 


GD40, GR384 


-6.02 


15319 


8.0 


0.587 


2.162 


5.866x10- 


13 


DBAZ"- 




HE 0413-3306 


04:15:20.63 


-32:59:10.0 


B 


= 15.9 




-4.81 


16640 


7.915 


0.543 


1.284 


3.088x10- 


12 


1 




HE 0414-0434 


04:16:52.65 


-04:27:23.6 


V 


= 15.7 




-5.20 


13 473 


8.0 


0.583 


3.379 


1.080x10- 


11 


new object^ 




WD 0503+ 147 


05:06:15.83 


+ 14:48:31.1 


V 


= 13.8 


KUV 05034+ 1445 


-4.91 


16246 


8.0 


0.589 


1.728 


3.289x10- 


12 


1.6 




WD 0853+ 163 


08:56:19.01 


+ 16:11:03.9 


mc = 


15.83 


PG 0853+164, LB 8827 


~ -3 


26128 


8.396 


0.853 








magnetic" 




WD 092 1+091 


09:23:55.27 


+08:57:17.5 


V = 


16.19 


PG 092 1+092, SDSS 


-4.81 


18713 


7.956 


0.571 


0.859 


6.753x10- 


4 






WD 0948+013 


09:51:02.21 


+01:04:32.1 


B = 


15.89 


PG 0948+013, SDSS 


-5.09 


16 939 


7.967 


0.572 


1.354 


1.061x10- 


12 


I 




WD1115+158 


11:18:22.77 


+ 15:33:33.3 


B = 


16.12 


PG 1115+158, DTLeo 


-2.44 


23 920 


7.523 


0.385 


0.198 


7.172x10- 


7 


V777 Hei- variable 




WD 11 34+073 


11:36:54.31 


+07:03:35.7 


B = 


16.32 


PG 1134+073 


-4.70 


18135 


8.315 


0.791 


2.304 


2.739x10- 


13 


DBAZ'-* 




WD 1149-133 


11:51:50.61 


-13:37:15.1 


V = 


16.29 


EC11492-1320, PG 


-4.12 


18618 


8.196 


0.714 


1.570 


4.807x10- 


13 






HE 1207-2349 


12:09:36.56 


-24:06:19.4 


B = 


15.79 




-5.24 


17 435 


7.935 


0.556 


1.107 


3.578x10- 


13 


new object 
















-5.26 


17 356 


7.921 


0.548 


1.092 


3.975x10- 


13 






EC 12438-1346 


12:46:30.42 


-14:02:40.8 


V = 


16.39 




-4.63 


16914 


7.933 


0.553 


1.253 


3.287x10- 


12 


1 




WD1311 + 129 


13:13:51.35 


+ 12:40:09.4 


mc = 


16.26 


PG 1311 + 129, LP 497-1 14, SDSS 


-2.55 


20224 


7.890 


0.539 


0.519 


4.878x10- 


14 


















-1.95 


22 403 


7.912 


0.555 


0.342 


2.573x10- 


15 






WD 1338-220 


13:41:17.29 


-22:19:43.7 


B 


= 16.2 


HE 1338-2204, EC 13385-2204 


-4.91 


13 694 


8.0 


0.584 


3.210 


1.925x10- 


11 


6 




HE 1349-2305 


13:52:44.26 


-23:20:06.5 


B = 


16.36 




-4.67 


18173 


8.133 


0.673 


1.494 


3.344x10- 


13 


DBAZ 




WD 1352+004 


13:55:32.38 


+00:11:24.1 


B = 


15.72 


PG 1352+004, SDSS 


-5.08 


14316 


8.0 


0.585 


2.752 


9.854x10- 


12 


DBAZ*" 




WD 1403-010 


14:06:19.96 


-01:19:32.6 


V 


= 15.9 


G 064-043, GR 272, SDSS 


-5.21 


15 784 


8.0 


0.588 


1.933 


2.536x10- 


12 


6 




HE 1409-1821 


14:11:48.67 


-18:35:05.6 


B 


= 15.6 




-5.14 


18514 


8.099 


0.654 


1.274 


5.681x10- 


14 


I 
















-5.32 


18321 


8.019 


0.606 


1.095 


5.678x10- 


4 






WD 1421-011 


14:24:29.21 


-01:22:16.7 


B = 


15.97 


PG 1421-011 


-4.08 


16788 


7.997 


0.589 


1.503 


1.220x10- 


11 


DBAZ'* 




WD 1444-096 


14:47:37.02 


-09:50:05.5 


V = 


14.98 


EC 14449-0937, PG1444-096 


-5.54 


16 724 


8.044 


0.617 


1.710 


4.066x10- 


13 


1 
















-5.58 


16622 


7.989 


0.584 


1.540 


4.712x10- 


13 






WD 1456+103 


14:58:32.74 


+ 10:08:17.6 


B = 


15.89 


PG 1456+103, CW Boo 


-3.42 


21533 


7.904 


0.549 


0.404 


2.443x10- 


16 


V777 Her variable 




WD 1542-275 


15:45:22.57 


-27:40:06.4 


V 


= 15.5 


LP 916-027 


-5.61 


10 826 


8.0 


0.577 


6.789 


7.422x10- 


12 


only He I 5876 A & Ha' - ' 


WD 1557+192 


15:59:21.08 


+ 19:04:09.0 


V 


= 15.4 


KUV 15571 + 1913 


-4.88 


18 690 


8.100 


0.654 


1.223 


7.099x10- 


4 


I 




WD 1709+230 


17:11:55.68 


+23:01:01.5 


V 


= 14.9 


GD205, BPM 92077 


-4.36 


18 527 


8.067 


0.635 


1.177 


3.303x10- 


13 


DBAZ 
















-4.41 


18 488 


8.108 


0.659 


1.309 


3.194x10- 


3 






WD 1917-077 


19:20:34.96 


-07:40:00.8 


V = 


12.28 


L923-021,LTT7658 


-5.16 


10195 


8.0 


0.576 


8.280 


2.069x10- 


11 


only He i 5876 A & 


Ha' 


WD 2 130-047 


21:33:34.86 


-04:32:24.9 


V 


= 14.5 


GD 233, L 1002-062 


-5.34 


17 562 


7.951 


0.565 


1.115 


2.266x10- 


13 


! 




WD2142-169 


21:45:27.60 


-16:43:08.0 


B 


= 14.6 


HK 22944-44 


> -2 


39138 


7.737 


0.505 








D0= 




WD 2 154-437 


21:58:01.55 


-43:28:00.8 


V = 


15.04 


L 427-060, BPM 44275 


-4.78 


16 734 


8.018 


0.601 


1.601 


2.472x10- 


12 


! 




WD 2229+ 139 


22:31:45.50 


+ 14:11:13.1 


mc - 


15.99 


PG 2229+139, GR 906 


-4.54 


15 423 


8.0 


0.587 


2.107 


1.607x10- 


11 


6 




HE 2237-0509 


22:39:40.75 


-04:54:17.2 


V 


= 14.0 




-4.80 


12 134 


8.0 


0.581 


4.731 


4.175x10- 


1 


only He i 5876 A & 


Ha' 


WD 2253-062 


22:55:47.43 


-06:00:51.4 


y - 


15.06 


GD 243, EG 230 


-4.24 


17 380 


7.981 


0.581 


1.256 


3.779x10- 


2 






WD 2316-173 


23:19:35.39 


-17:05:28.8 


V = 


14.04 


G 273-013, L0791-040 


-5.27 


10 868 


8.0 


0.577 


6.708 


1.625x10- 


11 


only He i 5876 A & 


Ho-' 


HE 2334-4127 


23:37:38.74 


-41:10:32.7 


B = 


15.61 


(LTT9631) 


-5.34 


18250 


8.029 


0.611 


1.140 


6.369x10- 


14 


new object 





' First reported as a DBA here; previously known as DB. 

^ Since WD 0853-1-163 is a magnetic DBA, the fit is uncertain because the Zeeman-splitted lines were fitted with a non-magnetic model. 

^ Temperatures and hydrogen abundances where determined from the equivalent widths of He i 5876 A and Ha, respectively. The surface gravity was assumed to be logg = 8.0. 

WD 1134-1-073 and WD 1421-011 show Can lines but are not mentioned among the SPY DBZ/DBAZ objects of Koester et al. ( I2005bb . 
' WD 2142-169 shows a He i-dominated spectrum with Balmer lines and Hen 4686 A. 

* Fit using the alternate set of models with van der Waals line broadening. The surface gravity was kept fixed at log g = 8.0 for these fits. 
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6520 6540 6560 6580 6600 

Fig. 7. The Ua line in the two SPY spectra of WD 0453-295. 
The lower spectrum shows two line cores. Taken from 
Napiwotzki et al. fi2005l) . 



in one of the spectra. Thus the DB component is indeed a DBA 
or perhaps even a DAB star, whose H abundance has to be con- 
siderable, since the equivalent widths of Ha of both the DA and 
its companion have to be comparable. 

4.3.3. HE 21 49-051 6 

The spectrum of HE 2149-0516 shows strong Balmer lines and 
weaker but strong He i lines. It is a new DAB star. A pure helium 
(or hydrogen) atmosphere fit does not give a result since the line 
strengths are too strongly affected by the other constituent of 
the atmosphere; only a pure hydrogen atmosphere fit that is lim- 
ited to the Ha line results in a reasonable temperature estimate 
of 29 000 K to 30 000 K, but with a very low fitted gravity of 
logg X 6.9. Alternatively, we find T^ff = 32 000 K if the gravity 
is kept fixed at logg - 8.0. These fits of the line are not good 
and the results are thus only very rough estimates. Nevertheless, 
the resulting temperature is similar to those that were derived for 
most other stars of this class, i.e., near the low-temperature edge 
of the DB gap. 

4.4. DB and DBA mass distributions 

Fig.|8]shows the mass distribution of the DB and DBA stars. The 
magnetic DBA WD 0853-H163, the cool DO, and the objects for 
which log g was kept fixed during the model atmosphere fit are 
not included in this plot. The overall shape of the distribution is 
similar to that which is found for the DA stars by, e.g., Liebert 
et al. ( |2005] |. Giovannini et al. (2005 ), or M adej e t al. ( [2004 t, 
as well as for the SPY DA stars (Voss et al. 120071 in prepara- 
tion). The mean mass of all stars is 0.596+0.072 Mq, and thus 
very similar to the average mass of the DA stars that has long 
been found to be about 0.55 - 0.6 Mq. The similarity is especially 
striking for the subsample of DBs that do not show any hydro- 
gen in their spectra; their mean mass is 0.584+0.059 Mq. The 
mean mass of the DBA subsample is only very slightly higher, 
0.607+0.082 Mq. Such a trend towards higher masses in DBA 
stars was already found by Beauchamp et al. (1996 ), although 
much more pronounced: The mean masses for their sample of 
13 DBA stars and 41 DB without spectroscopically visible hy- 
drogen are <M>= 0.642 Mq for the DBA and <M>= 0.567 Mq 



.8 

Mass [M3„„] 

Fig. 8. The mass distribution of the DB stars. The filled part of 
each bin represents the DBA stars. 

for the DB. They speculated that the presence of hydrogen could 
be favored by a larger stellar mass through an increased hydro- 
gen accretion from the interstellar medium, or through a thinner 
convection zone in which hydrogen dilution would be decreased. 

However, the differences that we find between our DB 
and DBA distributions are not significant. Furthermore, to test 
Beauchamp et al.'s speculation, we tried to identify any depen- 
dency of the hydrogen abundance, the hydrogen mass, or the hy- 
drogen accretion rate on the stellar mass. No such dependency is 
found, which adds confidence to the assumption that DBAs and 
DBs are not intrinsically different from each other. 

Beauchamp et al. also found that the distribution of the DB 
masses is narrower than that of the DA stars, only 5% of their 
DB stars are placed outside of a mass interval from 0.5 Mq to 
0.65 Mq. Here, a fraction of 25% of the DBs are found outside 
that interval, which is not too different from the 32% of the SPY 
DA stars that have masses outside this range. 

Finally, the DB mean mass of this analysis is much lower 
than the unexplainedly high value of 0.77 Mq that was found in 
paper I. The difference is probably due to the fact that paper I 
used an earlier, less sophisticatedly reduced version of the data. 
A major improvement was a better correction of the response 
function of the spectrograph in the blue spectra. The imperfect 
coiTection in the previous version was likely responsibe for a sig- 
nificant fraction of the systematic offset. The reason for the high 
paper I masses may be not entirely clear, but it is reassuring that 
we now find a mass distribution that agrees with that of the well 
studied 80% of all White Dwarfs, the DA. According to these 
results, there is no need to assume any distinction between the 
structures of white dwarfs with hydrogen- or helium-rich atmo- 
spheres, or, for that matter, between DB and DBA stars. 

4.5. Hydrogen abundance and total liydrogen mass in fhe 
DBA stars 

Fig. |9] shows that the observed hydrogen abundances are lower 
for cooler DBA stars than for the hottest objects. This decrease 
of the average abundance is not unexpected since the helium 
convection zone deepens with decreasing temperature, leading 
to increased dilution of any hydrogen mass originally present. 

The lowest observed Ha equivalent widths in our DBA sam- 
ple are close to 300 mA, and therefore we assume this value as 
the detection limit for hydrogen in our data. The detection limit 
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-7 I 1 1 1 1 1 1 1—1 

10000 12000 14000 16000 18000 20000 22000 24000 

Tetf[K] 

Fig. 9. The hydrogen abundances of the DBA stars as a function 
of the effective temperature. The tilted crosses show the objects 
for which the alternate model grids were used and for which the 
surface gravity was kept fixed during the fit. Asterisks show the 
stars for which temperatures and abundances were derived from 
the equivalent line widths of Hei 5876 A and Ha, respectively. 
The dotted line corresponds to the Ha equivalent width detection 
limit of 300 mA. 

in terms of abundance vs. temperature is shown by the dotted 
line in the figure. The observed abundances reach near this limit 
at temperatures around 1 8 000 K but are higher for the hottest 
and coolest objects. This reflects the average quality of our spec- 
tra which is somewhat lower for the hot and cool objects, and 
better at intermediate temperatures. 

A more instructive quantity is the hydrogen mass itself, 
which is shown in Fig. [10] as a function of temperature and in 
Fig. [TT]as a function of white dwarf age. Since we assumed a 
constant mass below approx. 16 000K, the qualitative similar- 
ity between the two curves is not too surprising. Both of them 
clearly show that the deduced H masses increase with the age 
of the WD. Before interpreting this in a too simplistic way one 
should note, however, the observational detection limit, which 
prevents any observation of H masses below 10"'^ at the 
cool end. On the other hand, the absence of larger H masses at 
the hot end is explained easily: A WD with H mass of 10 Mq 
at 22 000 K would have a H/He number ratio > 1, and thus look 
more hke a DAB than a DBA. 



-17 I 1 1 1 1 1 1 1 I 

10000 12000 14000 16000 18000 20000 22000 24000 

Fig. 10. The hydrogen masses plotted against the efifective tem- 
perature; the symbols and the dotted line are the same as in fig. 
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Fig. 11. The hydrogen masses plotted against log(Age); the sym- 
bols are again the same as in fig.|9] The distribution appears more 
scattered along the age axis than it is along the temperature axis 
in fig. [TOl This is due to the effect of stellar mass on the cool- 
ing rate. The dotted line depicts the hydrogen detection limit for 
log^ = 8.0. This 300mA-limit would be placed at higher ages 
for stars with \ogg > 8.0. 



5. Discussion and conclusions 

The detection of low hydrogen abundances in such a large num- 
ber of spectra significantly increases the number of DBA stars 
among the white dwarfs with helium-rich atmospheres. With a 
fraction of 55%, the majority of DBs in our sample are DBAs. 
The increase is not unexpected because it seemed always plau- 
sible that decreasing the detection limit would discover even 
weaker hydrogen lines. There are 6 new objects in our sample 
which have not been reported before in the literature, and 4 of 
these are DBAs. In addition to that, 14 stars of our sample cata- 
loged as DBs turn out to be DBAs. One of them is the V777 Her 
variable WD1115H-158. 

How do these observations fit with our understanding of DB 
evolution and the DB gap? The evolution of a 0.6 Mq DA in the 
DB gap and its transformation into a DB by convective mix- 
ing has been studied in great detail by MacDonald & Vennes 



( 119911 ). Depending on the convection description used, the tran- 
sition can occur between 25 000 and 32 000 K, if the hydrogen 
mass is 10"'^ Mq. For 10"''* Mq of H mass these values go down 
to 11 800 - 18 000K. Since we observe DBA white dwarfs up 
to 24 000 K in our sample, we assume x 10"'^ Mq as a reason- 
able starting value for those DA within the DB gap, which will 
turn into hot helium-rich stars around 30 000 K. Because of the 
shallow convection zone this leads to the high observed values 
of log H/He = -3 to -2; the apparently pure He DBs in this tem- 
perature range may have just slightly lower H masses. Because 
of the steep increase of the convective mass such small amounts 
of hydrogen would become unobservable after a short evolution 
to temperatures below 20 000 K. We still find hydrogen in these 
cooler stars, but now indicating much larger total H masses of 
> IO^'-'Mq. We can envisage two scenarios for the origin of 
this hydrogen: Mixing of a thicker original H layer, i.e. transfor- 
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Fig. 12. The hydrogen accretion rate, shown on a logarithmic 
scale against the effective temperature. The symbols are the same 
as in the previous figures. 




8.2 
log(Age/yr) 



to be normal DAs, and only careful observation and analysis re- 
veals their unusual composition. It is currently completely un- 
known how many more of these objects might be hidden among 
the "normal" DAs. Overall, our results thus strongly support the 
idea that small but varying amounts of hydrogen are accreted 
over their lifetimes by all DB white dwarfs. 

The recent study by Dufour et al. ( 2007 ) of DZ white dwarfs 
provides a very interesting continuation of our results towards 
lower temperatures. In the overlap region between 10000 and 
12 000 K, our results for the H/He ratios, total H masses, and 
necessary accretion rates are in very good agreement. With some 
precaution because of the selection effects, we can extend their 
findings of the decreasing H accretion rate with effective tem- 
perature much beyond their limit of 12000K up to 24 000 K. 
This may have some relevance for one of the remaining mys- 
teries, which is strongly re-emphasized by Dufour et al. (2007): 
In those cases, where metal traces are detected in addition to H 
(the DBZ, or DBZA), the accretion rate for H compared to e.g. 
Ca is orders of magnitude below that expected from solar abun- 
dances in the accreted matter If there is some kind of screening 
mechanism involved, which prevents or hinders the accretion of 
hydrogen, it seems to increase its efficiency in a rather smooth 
way towards high temperatures. This would argue against the 
propeller mechanism (Wesemael & Truran [T982| l, which is most 
often invoked as a possible explanation. The only other possi- 
bility we are aware of is the existence of a weak stellar wind in 
hot DBs, which dies down with decreasing temperature. In their 
wind model to explain the carbon abundance in DBs, Fontaine et 
al. (I2005I I assumed a linear relationship between wind strength 
and age, and that it dies at 20 000 K. Whether such a very weak 
wind (comparable in absolute value to our derived accretion 
rates) exists down to very much lower temperatures, and whether 
it can explain the difference in the Ca/H accretion rates remains 
to be studied. 

Acknowledgements. B.V. and D.K. acknowledge support by the Deutsche 
Forschungsgemeinschaft (DFG) under project numbers K0738/21-1, 
K0738/22-1, and K0738/23-1. 



Fig. 13. The hydrogen accretion rate, shown on a logarithmic 
scale against the age of the star The symbols are the same as in 
the previous figures. 



mation of a DA into a DB at much lower temperatures than the 
canonical 30 000 K, or interstellar accretion. 

In view of the MacDonald & Vennes (1991) calculations we 
can discard the first scenario. None of their different convection 
models allows mixing of 10"'^ Mq at a temperature higher than 
1 2 500 K, whereas such masses are observed already at 1 8 000 K. 
Assuming a constant accretion rate over the lifetime of the DB, 
the accretion rates necessary to explain the observations are ex- 
tremely low: Between 10"'^ and 10"^' Mq per year (see FigsfT2l 
and[T3Tl. In their important study of the diffusion/accretion sce- 
nario for DZ white dwarfs, Dupuis et al. (,1993 J estimate a rate of 
2.7x 10~^^ Mq for a typical WD in the thin, warm phase of the in- 
terstellar medium. This is based on the Eddington rate, using just 
the gravitational cross section of a star in classical mechanics. It 
is very hard to imagine a lower rate, unless some unknown mech- 
anism completely prevents accretion, but in any case such or 
even lower rates would not lead to observable hydrogen. On the 
other hand, a significantly higher rate - continuous or episodic - 
would lead to H/He ratios of 0.01 to 0.1 at the low temperature 
end, as observed in HS 0146+1847 (Koester et al. i2005al l and 
GD362 (Zuckerman et al. |2007D . These stars appear at first sight 
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Fig. A.l. The spectral regions around Ha for all DBA stars, continued in Figs. lA.2l and lA.3l The intensity is given in arbitrary 
units; the wavelengths are in A. Note that the wavelength range is not the same in all plots. 
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Fig. A.2. The spectral regions around Ha for all DBA stars, continued from Fig. lA.ll 
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Fig. A.3. The spectral regions around Ha for all DBA stars, continued from Figs. lA.ll and lA.2l 



